i/1 


-fiD-A142  772 
UNCLASSIFIED 


SOURCES  OF  ERROR  IN  OBJECTIVE  ANALVSI5CU)  NAVAL 
POSTGRADUATE  SCHOOL  NONTEREV  CA  R  FRANKE  HAV  84 
NPS52-84-0882 

F/G  12/1 


NL 


AD-A142  772 


NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY  CALIFORNIA  93943 


R.  H.  Shumaker 
Commodore,  U.  S.  Navy 
Superintendent 


D.  A.  Schrady 
Provost 


This  work  was  funded  by  the  Naval  Environmental  Prediction  Research 
Facility,  Monterey,  CA  under  Program  Element  61153N,  Project  (none), 
"Interpolation  of  Scattered  Meteorological  Data". 

Reproduction  of  all  or  part  of  this  report  is  authorized. 


Richard  Franke 
Associate  Professor  of 
Mathematics 


Reviewed  by: 


Released  by: 


L 


:T  M 


Gordon  E.  Latta,  Chairman 
Department  of  Mathematics 


Kneale  T.  Marshall 
Dean  of  Information  and 
Sciences 


SECURITY  CLASSIFICATION  of  THIS  PACE  (TW»«n  Palm  lmw<0 


REPORT  DOCUMENTATION  PAGE 


I.  REPORT  NUMBER  |2.  OOVT 

NPS-53-84-0003  /JO  *  t 

4.  title  (tnd  Submit)  I 

Sources  of  Error  in  Objective  Analysis 


a  Arc  READ  INSTRUCTIONS 

rwc  _ BEFORE  COMPLETING  FORM 

2.  OOVT  ACCESSION  NC  1.  RECIPIENT'S  CATALOG  NUMBER 


IT.  AUTHORf*; 


f  TYPE  OF  REPORT  A  PERIOD  COVERED 

Technical  Report 
Interim,  PY  1984 

•  PERFORMING  ORG.  REPORT  NUMBER 
I.  CONTRACT  OR  GRANT  NUMBER^ 


Richard  Franke 


I.  PERFORMING  ORGANIZATION  NAME  and  ADDRESS 

Naval  Postgraduate  School 
Monterey,  CA  93943 

it.  controlling  office  name  and  adores* 


to.  PROGRAM  ELEMENT.  PROJECT.  TASK 
AREA  A  WORK  UNIT  NUMBERS 

Program  element  61153N 


controlling  office  name  ano  address  *2.  report  date 

May  1984 

Naval  Environmental  Prediction  Research  Facility  .  . 

,,  n,n,  .  IS.  NUMBER  OF  PAGES 

Monterey,  CA  93943  .g 

MONITORING  AGENCY  NAME  A  AODRESSfff  5 IhSmS  fro*  Controlling  Oltlca)  IS.  SECURITY  CLASS,  (at  I him  rmpon) 


Naval  Air  Systems  Command  Ui 

Washington,  DC  20361 

"is* 

I*.  DISTRIBUTION  STATEMENT  (ol  thlm  Report)  — - 

Approved  for  public  release;  distribution  unlimited 


Unclassified 


<Sa.  OECL  ASSI FIC  ATION/  DOWNGRADII'  ' 
SCHEDULE 


jt7.  DISTRIBUTION  STATEMENT  (ot  the  abatract  an  (indln  Block  20,  li  ditto  rant  from  Report) 


is.  supplementary  notes 


I  IB.  KEY  WORDS  (Continue  on  reeerae  aide  It  nocooeory  end  Identify  by  block  nunbtrj 


Objective  analysis 
Optimum  interpolation 
Thin  plate  splines 


Barnes’  method 
Splines 

Statistical  interpolation 


20.  ABSTRACT  (Conilnu*  on  revere*  ml,m  If  naea.aar r  ■»<  IdenllfF  Ay  block  nuetber) 

The  error  in  objective  analysis  methods  that  are  based  on  corrections 
to  a  first  guess  field  is  considered.  An  expression  that  gives  a 
decomposition  of  the  error  into  three  independent  components  is  derived. 

To  test  the  magnitudes  of  the  contribution  of  each  component  a  series  of 
computer  simulations  was  conducted.  Grid-to-observation  point  interpolation 
schemes  considered  ranged  from  simple  piecewise  linear  functions  to  highly 
accurate  spline  functions.  The  observation-to-grid  interpolation  methods 


DO  1473 


EDITION  OF  I  NOV  SB  IS  OBSOLETE 

S/N  0102-  LF.  01 4- 6601 


UNCLASSIFIED _ 

SECURITY  CLASSIFICATION  OF  THIS  PAGE  fBPier.  DMA  Enter**; 


UNCLASSIFIED _ 

MCUMTV  CLASSIFICATION  OF  THIS  FAOS  0«a  toiin4 


'  considered  included  most  of  those  in  present  meteorological  use,  such 
as  optimum  interpolation  and  successive  corrections,  as  well  as  proposed 
schemes  such  as  thin  plate  splines,  and  several  variations  of  these 
schemes.  The  results  include  an  analysis  of  cost  versus  skill;  this 
information  is  summarized  in  plots  for  most  combinations.  The 
degradation  in  performance  due  to  inexact  parameter  specification  in 
statistical  observation-to-grid  interpolation  schemes  is  addressed. 

The  efficacy  of  the  mean  squared  error  estimates  in  this  situation 
is  also  explored. 


DTIC 


B 


S'N  0 T02-  LF-  014*  4601 


UNCLASSIFIED 


•scumrv  CLASSIFICATION  or  THIS  F AO«(T»li*n  Dim  Bilme 


1.0  Introduction 


The  purpose  of  this  investigation  was  to  test  the  relative 
importance  of  various  aspects  of  correcting  predicted  values  on  a 
grid  by  incorporating  information  from  observed  values  at  scat¬ 
tered  data  points.  Grid  and  observation  configurations  were 
patterned  after  those  routinely  available  over  North  America. 
Although  investigations  were  limited  to  the  univariate  objective 
analysis  methods,  I  believe  the  results  are  indicative  of  those 
that  would  be  achieved  in  the  more  general  case. 

Previous  investigations  on  the  error  contribution  of  various 
steps  in  the  objective  analysis  process  are  limited.  Koehler 
(1979)  separately  studied  the  errors  of  a  number  of  grid-to- 
observation  and  observat ion-to-gr id  interpolation  (approximation) 
routines.  He  noted  that  although  little  attention  is  typically 
paid  to  the  gr id-to-observation  interpolation  process  significant 
errors  may  be  caused  by  this  phase  of  objective  analysis.  While 
this  may  be  a  surprise  since  these  errors  are  usually  small 
compared  to  the  first-guess  errors  at  the  grid  points,  my  results 
further  demonstrated  that  the  contribution  to  overall  error  made 
by  the  gcid-to-observation  interpolation  process  should  not  bo 
ignored.  This  investigation  complements  recent  work  by  Seaman 
(1991)  regarding  the  accuracy  of  statistical  and  successive  cor¬ 
rection  schemes.  His  work  provides  expected  mean  squared  error 
estimates  for  these  schemes.  His  work  is  very  thorough  in  that 
it  provides  estimates  of  the  analysis  error  as  the  parameters  of 
the  first-guess  error  are  varied  while  holding  the  assumed  values 
constant,  and  vice-versa. 


In  Section  2  I  derive  a  generalized  expression  for  1 1  >  ^ 
overall  error  in  objective  analysis  which  leads  to  several 
observations.  In  Section  3  I  describe  the  simulation  method  and 
the  various  options  which  can  be  easily  handled.  In  Section  d  T 
present  the  results  of  the  simulations  and  discuss  their  implica¬ 
tions  with  regard  to  the  observations  made  in  Section  2. 

2.0  The  Form  of  the  Error  Term  in  Objective  Analysis 

My  setting  for  study  of  the  objective  analysis  process 
assumes  the  following: 

(i)  The  true  field  (function)  to  be  analyzed  is  H. 

(ii)  H  is  known  imperfectly  at  grid  points  through  a  "first- 
guess"  which  is  in  error  by  an  amount  to  be  denoted  by  g.  The 
error  is  a  normally  distributed  stationary  random  function  which 
has  a  certain  spatial  correlation  and  standard  deviation. 

(iii)  H  is  imperfectly  measured  at  observation  points  yield¬ 
ing  values  with  errors  o.  These  errors  are  independent  and 
normally  distributed  with  certain  standard  deviation. 

The  nature  of  the  errors  makes  it  only  possible  to  evaluate 
g  at  grid  points,  and  o  at  observation  points,  although  it  is 
sometimes  convenient  to  think  of  them  as  functions  rather  than  as 
sets  of  errors.  The  objective  analysis  process  consists  of  int¬ 
erpolation  of  the  first-guess  values  from  the  grid  to  th?  obser¬ 
vation  points  (by  a. linear  operator  designated  M)  followed  by 
interpolation  of  the  difference  between  the  observe!  and  first- 
guess  values  back  to  the  grid  point  (by  a  linear  operator  desig¬ 
nated  L)  as  a  correction  to  the  first-guess  values.  Denote  the 
error  in  the  entire  process  by  E,  then  the  final  approximation  is 
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H  +  E  =  Htg*L(H  +  o-M(Htg)) 

Let  m(H)  represent  the  error  in  the  approximation  of  ;i  by  M(H), 
then  M (H)  =  H  -  m(H).  Rearranging  and  simplifying  the  above, 
leads  to 

E=g+L(H+o-M  (H)  -  M (g) ) 

=  g  +  L(H  +  o  -  H  +  m  (H)  -  M(g)) 

=  g  +  L(o  +  n  (H) )  -  LM  (g)  , 

and  finally, 

E  =  L  (o)  +  Lm{H)  +  (g  -  LM  (g ) )  .  (1) 

Thus  the  error  is  made  up  of  three  parts.  The  term  L(o)  is 

dependent  on  the  'function'  o,  which  describes  instrumentation 
error  and  is  typically  not  controlable.  It  is  obviously  advan¬ 
tageous  to  have  o  small.  Since  the  values  of  o  are  assumed 
independent  and  random  it  is  desirable  for  L  to  be  a  smoothing 
operator.  The  second  part,  Lm(H)  is  within  our  control  and  the 
gr id-to-observation  point  interpolations  error  should  be  made 
small.  If  it  is,  then  interpolation  of  the  error  back  to  the 
grid  points  by  L  is  also  small,  assuming  this  smoothing  operator 
is  typical  and  does  not  magnify  the  error.  The  third  part  (g  - 
LM(g))  is  the  error  in  interpolation  of  the  first-guess  error  at 
the  grid  points  to  the  observation  locations  by  M ,  then  back  to 
the  grid  points  by  L.  While  it  is  possible  that  a  certain  sym¬ 
biosis  between  parts  could  occur,  the  goal  is  certainly  Cor  each 
interpolation  process  to  have  small  errors.  Ideally  the  operator 
L  should  be  a  left  inverse  of  the  operator  M,  although  this  is 
almost  certainly  impossible. 

Partitioning  the  error  in  this  way  shows,  for  example,  that 
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using  a  better  interpolation  process  from  the  grid  to  the  obser¬ 
vation  points  should  decrease  the  overall  analysis  error.  In 
certain  realizations,  of  course,  the  errors  may  tend  to  cancel. 
Since  the  three  terms  represent  uncorrelated  errors,  the  total 
error  variance  over  many  realizations  will  tend  to  be  the  sum  of 
the  individual  variances.  Thus,  decreasing  any  one  will  load  to 
statistically  smaller  error  variances. 

3.0  The  Computer  Simulation  Methods 

In  order  to  simulate  the  behavior  of  the  overall  error  under 
various  interpolation  processes  and  first-guess  error  assump¬ 
tions,  a  modular  computer  program  was  written  to  give  several 
options  for  the  different  processes.  This  made  it  possible  to 
test  a  large  number  of  combinations  of  methods  and  assumptions. 

In  general  terms,  the  process  simulated  consists  of  the 
following  steps: 

(i)  An  underlying  mathematically  defined  function 
describing  the  field  to  be  analyzed  is  evaluated  on  a 
grid  of  points. 

(ii)  "first-guess"  error  is  generated  from  normal  random 
deviates  with  a  pre-specif ied  standard  deviation  and  spatial 
correlation. 

(iii)  "Observed  values"  are  generated  by  evaluating  the 
field  to  be  analyzed  at  the  observation  points,  and  adding 
normally  distributed  uncorrelated  random  deviates  to  these 

values. 

(iv)  The  first-guess  values  at  the  observation  points  are 
obtained  by  one  of  several  interpolation  schemes. 


(v)  Based  on  the  difference  between  first-guess  mJ  obser¬ 
ved  values  at  the  observation  locations,  "corrected"  values  at 
tha  grid  points  are  obtained.  I  will  refer  to  the  corrected 
values  as  the  analysis  values. 

Most  of  the  simulations  were  done  with  two  different  grids 
and  observation  point  sets.  One  was  based  on  a  2.5°  grid  cover¬ 
ing  112.5°  W  to  32.5°  W  and  30°  N  to  5C°  N ,  with  117  =  13x9  grid 
points  and  36  observation  points  within  the  grid,  as  shown  in 
Figure  1.  The  other  was  based  on  a  5°  grid  covering  125°  W  to 
75°  W,  and  25°  N  to  60°  N,  with  S3  =  11x8  grid  points,  and  77 
observation  points  within  the  grid.  This  grid  and  the  observi- 
tion  locations  are  shown  in  Figure  2.  All  the  simulations  used 
were  univariate  analysis  methods  on  a  two-dimensional  field. 

This  simplification  was  necessacy  for  two  reasons.  The  first 
reason  is  that  the  generation  of  error  with  a  specified  spatial 
correlation  required  factorization  of  the  correlation  rruitrix  into 
the  product  of  a  lower  triangular  matrix  and  its  transpose.  The 
correlation  matrix  is  of  order  equal  to  the  number  of  grid  points, 
and  it  is  not  particularly  well  conditioned.  Incorporation  of 
multiple  levels,  a  large  grid,  or  correlated  multiple  variables 
was  therefore  not  possible.  The  other  reason  is  that  statistical 
results  required  that  numerous  realizations  be  simulated,  thereby 
limiting  the  time  available  to  do  the  computations. 

The  underlying  mathematically  defined  field  can  be  any  spec¬ 
ified  function.  The  height  field  test  function  used  is  the  on? 
given  by  Koehler  (1979)  and  also  described  in  Wahba  and 
Wendelberger  (1930).  The  input  parameters,  0Q  (the  location  of 


the  longitudinal  wave),  A9,  (amount  part  of  the  field  is  skewed 
log i tud i nal ly) ,  and  p  (the  pressure  for  the  height  field)  are 
easily  varied.  The  experiments  simulated  the  50'J  nb  height 
field,  using  fixed  or  randomly  varying  9Q  and  A e.  a  typical 
field  of  height  contours  generated  by  this  function  is  shown  in 
Figure  3.  First-guess  errors  had  a  nominal  standard  deviation, 
Tg,  of  30  m.  The  spatial  correlation  function  was  modeled  using 
exp((-d/cd))  ,  where  d  is  distance  (on  the  degree  grid),  and  c.j 
is  a  correlation  distance,  specified  as  10°.  I  have  used  degree 
measure  for  distance  rather  than  true  distance,  to  maintain  a 
rectangular  grid  of  first-guess  points.  This  resulted  in  a 
distortion  of  the  distance  varying  with  location.  The  observa¬ 
tion  errors  had  a  nominal  standard  deviation,  rQ,  of  1(5  m.  The 
observation  locations  approximately  correspond  to  the  North  Ame¬ 
rican  radiosonde  network  within  the  grids  being  used.  They  are 
shpwn,  along  with  the  grids,  in  Figures  1  and  2. 

The  output  consisted  of  mean,  root-mean-square,  and  maximum 
errors  over  each  data  set  (first-guess  at  grid  points,  first- 
guess  at  observation  locations,  observation  at  observation  loca¬ 
tions,  and  analysis  values  at  grid  points)  for  each  realization. 
The  first  and  third  of  these  mainly  served  as  a  check  on  the 
psuedo-random  number  generator  (IMSL  subroutines  GGNS  >1  a  net 
GGN*IL).  The  output  also  gave  summaries  of  the  same  errors  over 
all  realizations  as  well  as  the  mean  and  standard  deviation  of 
the  root-mean-square  errors  over  the  realizations.  Interpolation 


processes  are  sometimes  ill-behaved  around  boundaries.  Since  in 
the  global  problem  this  can  be  avoided,  the  effects  were  mini¬ 
mized  here  by  tabulating  error  only  over  the  interior  grid 


points.  Thus  the  results  are  over  77  grid  points  on  the  ’T0 
grid  and  54  grid  points  on  the  5°  grid.  The  options  simulate.: 
for  each  step  are  described  below. 

a.  Grid-to-observation  point  interpolation 

First-guess  values  at  the  observation  points  are  obtained  by 
interpolation  from  the  first-guess  grid  values.  1  compared  four 
schemes.  Others  could  be  easily  included,  however  my  results 
indicate  it  will  probably  not  be  fruitful  to  do  so.  The  methods 
I  have  used  are: 

(i)  Piecewise  bilinear  interpolation.  As  with  any  piece¬ 

wise  defined  method,  one  must  first  determine  the  rectangle  in 
which  the  evaluation  point  lies.  Then,  the  evaluation  is  most 
easily  seen  as  translation  to  the  square  [0,1]  ,  followed  by  3 

one  dimensional  interpolations.  This  requires  8  operations, 
where  an  operation  is  defined  as  a  multiplication  or  division 
followed  by  an  addition  or  subtraction.  Practically,  the  evalua¬ 
tion  can  be  accomplished  in  5  operations  (and  a  couple  of  extra 
additions/subtractions).  In  my  cost  analysis  I  have  used  8 
operations;  this  cost  is  very  low  compared  to  that  of  other 
necessary  calculations. 

(ii)  Bicubic  spline  interpolation.  T  used  the  1M8L  subrou¬ 
tines  IBCCCU  and  IBCEVL.  Preprocessing  for  the  spline  coeffi¬ 
cients  on  a  NqxN^  grid  requires  1  2NqxN^  +  27:'1(j)  +  51  Ng-10  operations. 
Evaluation  requires  2  operations  to  translate  to  [0,1]  *  and  5 
cubic  interpolations  at  9  operations  each  for  a  total  of  47 
operations.  The  preprocessing  operations  involve  solution  of 
tridiagonal  systems  of  equations  which  ate  amenable  to  vector izn- 
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summarize  the  results  in  Table  l,  along  with  son1  r  ?p  r vs  -mi  ►  a  r  i  / 
numbers  that  arise  from  my  simulations.  I  have  describe-.!  some 
schemes  as  local,  implying  that  others  are  global.  In  the  con¬ 
text  of  global  objective  analysis,  all  the  schemes  1  consider  are 
local;  the  schemes  which  are  global  for  my  simulation  are  less 
local  than  the  ones  I  refer  to  as  local. 

(i)  Optimum  interpolation  (01).  This  scheme  was  i  nt  rodu-'-'J 
to  the  meteorological  literature  by  Gandin  (19e3)  .and  has  re¬ 
ceived  widespread  attention  in  recent  years,  e.g.  see  Bergman 
(1979)  and  Lorenc  (1981).  The  method  in  its  proper  form  requires 
that  the  spatial  covariance  function  of  the  first-guess  errors 
and  the  standard  deviation  of  the  observation  error  be  known. 
Since  these  are  known  for  this  simulation,  I  have  used  their 
properties.  I  have  implemented  the  scheme  as  described  in  Franke 
and  Gordon  (1983),  viewing  the  approximation  as  a  linear  combina¬ 
tion  of  the  covariance  function'-  associated  with  the  observation 
points.  Thus  we  have 
No 

AH(P)  =  E  akC(P,Pk)  , 
k  =  1 

where  C(P,Q)  is  as  noted  above.  The  ak  are  determined  from  the 
system  of  equations 

(  C  (P  i ,  P  j )  +  6ijr2) 

where  AHj  is  the  difference  between  the  first-guess  and  observed 
values  at  the  ith  observation  point,  rQ  is  the  standard  deviation 
of  the  observation  error,  and  is  the  Kronecker  delta. 

The  cost  of  (01)  consists  of  a  preprocessing  phase  that 


includes  the  generation  and  solution  of  th.*  system  of  equations, 
followed  by  evaluation  of  the  analysis  at  the  grid  points.  Fot 
Nq  observations,  preprocessing  is  at  a  cost  of  N  0  (  n'  Q  f  1 )  /  ?  func- 
tion  evaluations  to  generate  the  coefficient  matrix  and  (N~tSN0- 
N0)/6  operations  plus  NQ  square  roots  to  perform  Cholesxy  decom¬ 
position  and  solution  of  the  system  of  equations  for  the  ak. 
Evaluation  costs  NQ  covariance  function  evaluations  and  d0  opera 
tions  to  form  the  linear  combination  representing  the  value  of 
the  correction  at  each  grid  point. 

(ii)  Local  optimum  interpolation.  In  my  version  of  this 
scheme,  nominally  only  points  within  the  surrounding  13°  square 
are  used;  if  fewer  than  4  observations  are  available,  the  squar 
is  expanded  to  15°  and  so  on,  by  2.5°  increments  in  each  direc¬ 
tion  until  at  least  4  observations  are  available.  The  costs  of 
the  search  were  not  assessed.  For  each  grid  value  correction,  a 
system  of  equations  must  be  formed  and  solved,  and  the  correspon 
ding  correction  computed.  With  n  observations  being  used  the 
expressions  given  for  01  above  apply  with  n  replacing  mq.  This 
process  was  performed  for  each  grid  point,  making  the  total  cost 
the  sum  of  these  costs  over  all  grid  points. 

(iii)  Global  Barnes'  method.  This  type  of  scheme  is  des¬ 
cribed  by  Barnes  (1971)  and  others.  My  scheme  used  the  known 
correlation  functions  as  the  weights  for  the  first  pass.  Thus, 

the  approximation  is 

N  N 

o  o 

AH%1(P)  *  Z  Vvk(P)AHk/  I  Wk(P)  , 
k= 1  X= 1 

where  Wk  *  exp(  (- ( I  I  P-Pk  I  I /c^)  2) ,  and  A  Mk  is  as  before.  For  the 


second  pass  the  correction  has  the  same  form,  but  AHj,  is  r  ■  •  p  1  ■ .» :•  e  d 
by  the  difference  between  the  corrected  first-guess  and 

the  observations.  The  quantity  c^  is  replaced  by  c^/3^2  for  the 
second  pass.  The  total  correction  at  the  grid  points  is  then  the 
sum  of  the  two  corrections.  For  each  grid  point  the  cost  of  this 
method  is  NQ  weight  function  evaluations  per  pass  and  N Q 1  opera¬ 
tions  per  pass.  In  addition  a  separate  interpolation  from  the 
grid  points  to  the  observation  points  is  required  before  the 
second  pass.  This  type  of  scheme  has  been  defined  and  studied  in 
a  different  context,  without  a  change  of  weight  functions  between 
iterations,  by  Foley  and  Nielson  (193d). 

(iv)  Local  Barnes'  method.  The  same  localization  process 
as  used  for  the  local  01  scheme  (ii)  was  used  here.  As  for  the 
global  version,  two  passes  were  used.  Hence  the  cost  for  an 
evaluation  at  a  grid  point,  with  n  neighboring  points  is  the  same 
expression  as  in  the  global  scheme,  but  with  n  replacing  NQ.  In 
addition,  there  was  the  search  cost  to  determine  the  nearby 
observations,  which  was  not  assessed.  Costs  of  an  interpolation 
from  the  grid  points  to  the  observation  points  between  passes  was 
included. 

(v)  Statistical  interpolation  (c^  =  14°).  In  practical 
applications  of  01  the  error  correlations  and  standard  deviations 
cannot  be  modeled  precisely.  This  has  lead  to  the  use  of  the 
name  "statistical  interpolation".  Computationally  the  method  is 
identical  to  the  01  scheme  (i).  Here  the  only  difference  is  the 
substitution  of  an  inexact  correlation  distance,  c(J  -  14.  The 
algorithm  and  costs  are  identical. 

(vi)  Statistical  interpolation  (c^  =  7°).  Agiin  this  is 


identical  to  (i)  except  that  the  inexact  value  substitute'5  for 
is  7. 

(vii)  Statistical  interpolation  (damped  cosine  correlation 
function).  Once  more  this  scheme  is  computationally  identical  to 
(i)  except  that  the  correlation  function  used  is  of  the  form 
exp((-(|  IP-QI  l/cd)2)cos((|  IP-QI  l/cd) (  / 2 ) ) .  I  used  the  value  cd 

-  in. 

(viii)  Thin  plate  splines.  This  method  is  described  by 
Wahba  and  Wendlelberger  (1980)  and  others.  The  approximating 
function  used  by  the  scheme  is 


H  ( P )  =  Z  AkB(P,Pk)  +  ad  +  b<*  +  c  , 
k  =  l 

where  the  basis  function  B(P,Q)  =  II P-Q I  I 2log I  I P-Q ||  .  The  Ak 
and  a,  b,  and  c,  are  obtained  by  solving  the  system  of  equations 


N° 

Z  Ak(B(Pi,Pj)+XNkro  t  j)  +  aai+b(*  t+c  =  AHj  ,i»l,...,N( 

j-1 

^o 

Z  AjSj  =  0 

j-1 

No 

Z  Aj<kj  =  0 

j-1 

No 

Z  Aj  =  O  . 

j*l 


In  the  above,  X  is  a  smoothing  parameter.  The  smoothing  parame¬ 
ter  was  chosen  on  the  basis  of  a  few  trials  with  no  attempt  to 
optimize  its  choice  for  a  particular  data  set,  as  can  be  done. 
Wendelberger  (1931)  describes  a  program  that  will  automatically 


choose  X  (and  m  as  wel),  see  next  method),  but  I  have  not.  tested 


but:  nor  pos  i  t  i  ■■■'» 


it  yet.  This  system  of  equations  is  symmetric, 
definite.  I  have  used  standard  L-U  decomposition  routines  to 
solve  the  system.  Methods  for  symmetric  indefinite  systems  us” 
about  half  as  many  operations,  however  I  observed  greater  numeri¬ 
cal  stability  using  the  general  decomposition  process.  There  ere 
N0(N0+l)/2  basis  function  evaluations,  and  solution  of  the  system 
of  equations  requires  (MQ+3)  (Nq+^Nq+3) /3  +  (NQ  +  3)2  operations. 
Unlike  symmetric  positive  definite  systems,  solution  of  these 
equations  requires  searching  for  a  pivot  and  pivoting.  Evalua¬ 
tion  at  each  grid  point  then  requires  NQ  basis  function  evalua¬ 
tions  and  Nq+2  operations  to  form  the  sum. 

(ix)  Laplacian  smoothing  spline  (m=3).  This  scheme  is  also 
described  by  Wahba  and  Wendelberger  (198?.),  and  is  one  of  those 
available  in  the  program  by  Wendelberger  (1981).  The  thin  plate 
spline  method  is  a  member  of  this  family  (with  m  =  2),  but  also  has 
the  "thin  plate"  interpretation.  The  reason  for  inclusion  of 
this  method  is  that  the  results  of  Wahba  and  Wendelberger  indi¬ 
cate  that  pressure  height  fields  are  better  approximated  using 
values  of  m  =  3  or  4.  I  will  not  describe  the  method  fully.  It 
requires  evaluation  of  N0(NQ+l)/2  basis  functions  and  3NQ  multi¬ 
plications  to  set  up  the  system  of  NQ  +  5  equations  to  be  solved. 
Then  NQ+5  operations  would  be  required  for  evaluation  at  each 
grid  point,  along  with  tne  evaluation  of  basis  functions. 

(x)  Franke/Gordon.  This  scheme  was  suggested  by  Franks  and 
Gordon  (1983)  as  one  which  is  an  explicit  scheme,  similar  to 
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(in  the  notation  of  that  report)  were  performed.  The  :’ost  in 
operations  is  2NQ(N0+1)  plus  3NQ  for  each  grid  point.  The  number 
of  weight  function  evaluations  is  2N£  plus  3MQ  for  each  grid 

point. 

(xi)  Pseudo-Barnes’  method.  This  method  was  described  in 
Franke  and  Gordon  (1983)  and  was  at  that  time  mistaken  for  Barnes' 
method.  It  differs  in  that  the  error  at  the  second  iteration  is 
Barnes'  approximation  evaluated  at  the  observation  point  minus 
the  first-guess  error,  rather  than  the  the  corrected  first-guess 
at  the  grid  point  interpolated  to  the  observation  point,  minus  the 
observed  value.  The  cost  of  this  algorithm  is  evaluation  of 
weight  functions  plus  2NQ  for  each  grid  point.  It  requires 
N0(N0+1)  operations,  plus  2(NQ+1)  for  each  grid  point. 

(xii)  Local  pseudo-Barnes'  method.  This  is  a  local  version 
of  (xi),  using  the  same  "nearby"  observation  points  as  (ii)  and 
(iv).  A  grid  point  with  n  nearby  observation  points  requires 
evaluation  of  n£+2n  basis  functions  and  n‘+3n+2  operations. 

4.(1.  Results 

The  simulation  program  described  in  the  previous  section  was 
run  for  a  substantial  number  of  different  options.  Each  run 
consisted  of  100  realizations  of  a  test  field  each  containing 
associated  first-guess  and  observation  errors.  Table  2  gives  the 
assumed  parameter  values  for  the  various  cases.  Not  all  combina¬ 
tions  of  grid-to-observation  point  and  observation-to-gr id  point 
interpolation  schemes  were  used  in  every  case.  The  tables  detail 
the  complete  results  and  the  entries  indicate  which  combinations 
were  computed.  Each  combination  in  a  given  table  (3-14) 
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corresponds  to  the  same  set  of  real  izations,  but  different  tables 
depend  on  different  realizations. 

This  investigation  was  designed  to  determine  the  influence 
of  the  gr id-to-observat ion  point  interpolation  scheme.  This 
influence  is  seen  by  noting  changes  in  error  for  a  particular 
observation-to-gr id  point  interpolation  scheme  as  the  grid-to- 
observation  point  interpolation  scheme  is  varied.  The  rows  of 
Tables  3-14  give  this  information.  The  bicubic  spline 
interpolation  produced  significant  improvement  over  piecewise 
bilinear  interpolation.  This  verifies  the  smaller  magnitude  of 
the  term  Lm(H)  in  the  error  expression  given  by  (1)  for  the 
spline  method.  For  2.5°  grids  the  errors  were  no  smaller  for 
spline  interpolation  than  for  piecewise  bicubic  or  Bessel  bicubit' 
interpolation.  Evidently  the  grid  spacing  was  small  enough  (for 
the  test  function  used)  that  the  interpolation  error  was  not 
significant.  Spline  interpolation  did  show  an  improvement  over 
piecewise  bicubic  and  Bessel  bicubic  interpolation  on  the  5° 
grid.  Spline  interpolation  and  the  cubic  interpolation  methods 
showed  even  greater  improvement  over  piecewise  linear  interpola¬ 
tion  on  the  5°  grid  than  on  the  2.5°  grid.  Interestingly  the 
first-guess  errors  at  the  observation  points  had  greater  rms 
values  for  cubic  interpolation  than  they  did  for  linear  interpo¬ 
lation.  This  occurs  because  linear  interpolation  inherently  has 
greater  smoothing. 

’'lost  of  the  useful  information  given  in  Tables  3-11  can  be 
more  easily  obtained  from  plots  of  the  salient  values.  Figures 
4-8  give  plots  of  skill  vs.  cost  of  the  algorithm  in  thousands  of 


operations  per  analysis. 


Her 


"skill"  is  defined  to  bn  1 


rmsa/r0,  where  rmsa  is  the  rms  error  in  the  analysis  values.  The 
skill  with  respect  to  bilinear  and  bicubic  interpolation  are  each 
indicated,  connected  with  a  straight  line  to  delineate  the  extent 
between  the  two.  The  results  for  only  one  of  the  statistical 
schemes,  (vi),  has  been  plotted  since  the  others  were  nearly 
identical.  For  these  purposes  I  counted  an  evaluation  of  a  basis, 
weight,  square  root,  or  covariance  function  as  10  operations. 

The  plots  reveal  that  the  statistical  schemes,  local  01,  and  thin 
plate  splines  all  had  close  to  the  same  accuracy  and  all  were 
slightly  less  accurate  than  01.  The  Barnes'  schemes,  the 
Franke/Gordon  scheme,  and  Laplacian  smoothing  splines  were  least 
accurate.  The  poor  performance  of  the  Laplacian  smoothing 
splines  here,  in  contrast  to  the  better  performance  obtained  by 
Wahba  and  Wendelberger  (1930)  is  probably  due  to  the  scheme  being 
applied  to  the  first-guess  error  function  rather  than  to  the 
underlying  true  height  field.  The  degradation  in  the  performance 
of  the  less  than  optimal  statistical  schemes  is  perhaps  less 
drastic  than  one  might  expect.  It  does  appear  that  it  was  better 
to  underestimate  the  correlation  distance  than  to  overestimate 
it. 

Figure  9  shows  plots  of  the  rms  errors  in  the  analysis 
values  as  a  function  of  first-guess  errors.  The  improvement  in 
the  Barnes'  scheme  as  the  first-guess  errors  decrease  was  rapid. 
The  scheme  gave  results  nearly  as  good  as  01,  the  statistical 
schemes,  and  thin  plate  splines.  This  occurred  because  the 
principal  problem  became  smoothing  observation  errors  as  the 
first-guess  errors  tended  to  zero.  Figure  1G  shows  plots  of  the 

ir> 


rms  errors  in  the  analysis  values  as  a  function  of  observation 


errors.  As  observation  errors  go  to  zero  the  importance  of 
modelling  the  first-guess  error  was  more  important  than 
smoothing.  Thus  01,  the  statistical  schemes,  and  thin  plate 
splines  improved  the  most,  while  both  Barnes'  schemes  improved 
little.  Figures  11-13  show  the  rms  errors  in  the  analysis  values 
when  incorrect  variances  wore  specified  for  thn  interpolation 
routines.  Methods  not  using  these  values  were  naturally  unaf¬ 
fected  so  that  changes  in  the  rms  errors  in  the  analysis  values 
for  these  methods  only  reflect  the  variability  of  the  (different) 
realizations  used  in  the  various  cases.  The  plots  show  that  the 
use  of  incorrect  values  for  the  first-guess  and  observation  error 
variances  did  not  drastically  affect  the  accuracy  of  the  statist¬ 
ical  methods.  The  interested  reader  is  referred  to  Seaman  (198  5) 
for  more  extensive  tests  of  the  effects  of  incorrect  parameter 
specification  on  the  performance  of  statistical  interpolation 
methods. 

One  of  the  attractive  features  of  the  statistical  schemes  is 
that  they  afford  a  calculation  for  the  estimated  mean  squared 
error.  These  estimates  do  not  depend  on  any  particular  realiza¬ 
tion,  so  they  were  not  incorporated  into  the  process.  However,  I 
did  compute  them  as  a  side  calculation  for  my  grids  and  observa¬ 
tion  points.  The  results  of  these  calculations  are  tabulated  for 
the  2.5°  grid,  along  with  tnc  empirical  rms  errors  obtained 
during  the  simulations.  Table  15  shows  that  the  estimates  given 
by  01  were  quite  good;  the  estimated  and  empirical  errors  varied 
only  a  few  percent.  They  also  were  accurate  for  local  01,  as 
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they  should  be.  On  the  other  hand,  the  slight  degradation  in 
performance  of  statistical  methods  when  incorrect  correlations  oi 
variances  were  specified  did  not  carry  over  to  the  error  esti¬ 
mates.  In  fact  the  schemes  that  have  their  performance  degraded 
the  most  (in  this  case,  using  too  long  a  correlation  distance) 
showed  a  decrease  in  the  estimated  error  variance.  Conversely, 
shortening  the  correlation  distance  in  the  statistical  method 
increased  the  error  estimate  as  well  as  the  empirical  error 
obtained,  although  the  empirical  error  is  underestimated.  This 
indicates  that  one  must  not  put  too  much  faith  in  the  error 
estimates  when  the  actual  covariance  structure  is  not  known,  as 
in  practice.  It  appears  one  could  obtain  just  about  any  error 
estimate  wished  simply  by  specifying  unrealistic  parameters  for 
the  covariance  structure. 

The  principal  results  of  this  study  were  as  follows.  The 
decomposition  of  the  error  into  independent  components  in  (1) 
identified  possible  ways  to  decrease  the  analysis  error.  This 
lead  to  the  results  showing  the  contribution  of  the  grid-to- 
observation  interpolation  process,  the  necessity  of  smoothing  in 
the  observation-to-gr id  interpolation  process,  along  with  accuracy 
The  simulations  provided  confirmation  of  the  above  and  yielded 
information  concerning  the  sensitivity  of  statistical  interpola¬ 
tion  schemes  to  inexact  parameter  specification.  The  error  esti¬ 
mates  provided  by  statistical  schemes  were  shown  to  be  sensitive 
to  inexact  parameter  specification. 
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10 
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0 
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30 

10 

100 

13.775 

13x9,2.5° 

20 

10 

100 

0 

13x9,2.5° 

30 

5 

100 

a 

13x9,2.5° 

30 

10 

random3 

random*3 
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30 

10 
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( r 
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10 
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(r 

30 

10 
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(r 

30 

5 
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random*3 

13x9,2.5° 

(r 

5 

10 

random3 

random*3 

13x9,2.5° 

30 

0 

random3 

random*3 

13x9,2.5° 

30 

10 

random3 

random*3 

11x8,5° 

Table  2 


0Q  uniformly  distributed  in  (-82 . 5° , 1 12 . 5°) 

A0  uniformly  distributed  in  (-15°, 15°) 

Tne  statistical  interpolation  routines  were  given 
incorrect  variances,  as  indicated 


W 


F  =  50C.  lheta=  100,  Eelth  =  0 

Entries:  RMSE 

Analysis 

rg  =  3  C ,  ro  =  10 

Kunber  cf  realizations  = 

100 

Mean 

RMSE  (StD®v) 

13X9  grid  of  2.5  degrees. 

36  observation 

points 

Grid-to-cbs:  FW  linear 

Bicub  Spl 

PW  3icub 

Bsl  Bicub 

Obs-to-grid 

Opt  Inteip 

6.64 

6  .09 

6.09 

6.09 

(Cd  =  10) 

6.53  (1.  18) 

5 .97  (1.20) 

5.98  (1.  19) 

5.98 (  1. 19) 

Local  OI 

7.09 

6.53 

6.54 

6.55 

(Cd  =  10) 

6.99  (1.19) 

6.42  (1.  19) 

6.44  (1.  19) 

6.44  (  1.  18) 

Barnes* 

9.27 

8.87 

8.87 

8.88 

2- Pass 

9.08  (1.03) 

8.68  (1.82) 

8.68  (1.82) 

8.69  (1.82) 

Barnes' 

8.42 

7.95 

7.96 

7.96 

(Lccal) 

8.27  (1.  57) 

7.79  (1.56) 

7.80  (1.56) 

7.81  (1.56) 

Stat  Interp 

7.28 

6  .78 

6.79 

6.79 

(Cd  *  14) 

7.23(  1.22) 

6.66  (1.27) 

6.67  (1.26) 

6.68  (1.26) 

Stat  Interp 

7.34 

6.87 

6.87 

6.87 

(Cd  =  7) 

7.23(1.25) 

6  .75  (1.26) 

6.76  (1.25) 

6.76 (1.25) 

Stat  InterF 

7.37 

6.91 

6.91 

6.91 

(DmFd  Cos) 

7.26  (1.28) 

6  .79  (1.28) 

6.79  (1.28) 

6.79  (1.28) 

Thin  PI  Spl 

7.12 

6.59 

6.60 

6.60 

(m  *  2) 

7.00  (1.30) 

6  .45  (1.33) 

6.46  (1.32) 

6.47  (1.32) 

Lapl  Sot  Spl 

10.54 

10.25 

10.25 

10.25 

(m  =  3) 

10.40  (  1.  73) 

10.10  (1.73) 

10.10  (1.73) 

10.11  (1.73) 

Frnke/Gr dn 

12.02 

11  .75 

11.76 

11.76 

(3  Pass) 

11.72(2.  65) 

11 .45(2.65) 

11  .45  (2.65) 

11.45  (2.65) 

PseudoBarnes* 

9.28 

8.87 

8.87 

8.88 

(2  Fass) 

9.10  (1.  83) 

8  .68  (1.82) 

8.68  (1.82) 

8.69  (1.82) 

PseudoBarnes' 

8.20 

7.70 

7.71 

7.72 

(Lccal) 

8.06  (1.51 ) 

7.55  (1.50) 

7.57  (1. 50) 

7.57(1.  50) 

TABLE  3 

RMSE  Analysis 
Mean  RMSE (StDev) 


E  =  500,  Theta  =  100,  Delth  *  13.775  Entries: 

rg  =  3C,  to  =  10 

Humber  of  realizations  *  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


Gr id-to-obs : 

Obs-to-grid 

EH  Linear 

Bicub  Spl 

PH  Bicub 

Bsl  Bicub 

Opt  Inteip 

6.88 

6.19 

6.20 

6.21 

(Cd  =  10) 

6.72(1.49) 

6 .05  (1.34) 

6.06  (1.34) 

6.06  ( 1.  34) 

Local  OI 

7.40 

6.76 

6.77 

(Cd  =  10) 

7.24  (1.51) 

6  .62  (1.37) 

6.63  (1.37) 

Barres' 

9.73 

9.28 

9.29 

2- Eass 

9.53(1.96) 

9.09  (1.88) 

9.10  (1.88) 

Barnes  * 

8.29 

7  .79 

7.77 

(Lccal) 

8.  16  (  1.  48) 

7.67  (1.36) 

7.66  (1.36) 

Stat  Interp 

7.71 

7.08 

7.09 

(Cd  =  14) 

7.56  (1.50) 

6.95  (1.35) 

6.96  (1.35) 

Stat  Interp 

7.54 

6.96 

6.96 

(Cd  =  7) 

7.39  (1.50) 

6  .82  (1.36) 

6.83  (1.36) 

Stat  Interp 

7.37 

6  .91 

6.91 

6  .91 

(D*pd  Cos) 

7.26  (1.28) 

6  .79  (1.28) 

6.79  (1.28) 

6.79  (1.28) 

Thir.  PI  Spl 

7.45 

6.80 

6.81 

(m  =  2) 

7.29  (  1.  51) 

6.66  (1.37) 

6.67  (1.37) 

Lapl  Sb  Spl 
(B  =  3) 


Frr.ke/Grdn 
(3  Eass) 

PseudoBarnes*  9.75  9.28 

(2  Eass)  9.55  (1.  97)  9.09  (1.88) 

PseudoBarnes* 

(Local) 


TABLE  4 


p  =  500,  Theta  =  100,  Celth  *  0  Entries:  RUSE  Analysis 

ig  =  20,  tc  =  10  Wean  RMSE(StDev) 

Number  of  realizations  =  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


Grid-to-cts : 

obs-to-grid 

PM  Linear 

Bicub  Spl 

PH  Bicub 

Opt  Inteip 

6.23 

5.75 

5.76 

(Cd  *  10) 

6.10  (1.  28) 

5.62  (1.22) 

5.63  (1.22) 

Local  OI 

6.54 

6.10 

6.10 

(Cd  =  10) 

6.4  1  (1.28) 

5.97  (1.22) 

5.98  (1.22) 

Barnes' 

7.18 

6.85 

6.86 

2- Pass 

7.03  (1.47) 

6.71  (1.47) 

6.71  (1.41) 

Barnes ' 

7.19 

6  .77 

6.76 

(Lccal) 

7.08(  1.22) 

6.68(1.11) 

6.67  (1.10) 

Stat  Interp 

6.70 

6.30 

6.31 

(Cd  =  14) 

6.58  (  1.27) 

6.19(1.19) 

6.19  (1.  19) 

Stat  Interp 

6.71 

6.24 

6.25 

(Cd  =  7) 

6.58  (1.  30) 

6.12(1.25) 

6.12(1.24) 

Stat  Interp 

6.78 

6.31 

6.32 

(Ompd  Cos) 

6.65  (1.33) 

6.18  (1.28) 

6.18  (1.28) 

Thir  PI  Spl 

6.66 

6.20 

6.20 

(n  *  2) 

6.52  (1.34) 

6 .06  (1.29) 

6.07  (1.29) 

Lapl  Sm  Spl 

11.05 

10.71 

10.71 

(m  =  3) 

1C. 84  (2. 16) 

10.50  (2.  12) 

10.50  (2.12) 

Frnke/Grdn 

8.72 

8.56 

8.56 

(3  Pass) 

8.54  (1.78) 

8.39  (1.70) 

8.39  (1.  70) 

PseudoBarnes' 

7.19 

6.85 

6.86 

(2  Pass) 

7.03(1.47) 

6.85  (1.41) 

6.71  (1.41) 

Bsl  Bicub 


PseudoBarnes'  6.78 
(Local)  6.66  (  1.28) 


6 .37 

6.25(1.24) 


6.38 

6.26  (1.23) 


f  =  500,  Theta  =  100,  Celth  =  0  Entries: 

rg  *  3C,  to  =  5 

Dasher  of  realizations  *  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


BMSE  Analysis 
Mean  BMSE(SxDev) 


Grid-to-cfcs:  EM  Linear  Bicufc  Spl  PM  Bicub 


Obs-to-grid 


Opt  Interp 
(Cd  =  10) 

4.57 

4.50  (0.83) 

3.76 

3.70  (0.69) 

3.77 

3.70  (0.69) 

Local  01 
(Cd  =  10) 

5.03 

4.95  (0.91 ) 

4.26 

4 .19  (0.77) 

a  27 

4. *20  (0.7  7) 

Barnes' 

2~Eass 

8.88 

8.69  (1.85) 

8.48 

8  .29  (1.76) 

8.49 

8.30  (1.77) 

Barnes ' 

(Local) 

6.39 

6.30  (  1.  10) 

5.90 

5.79(1.11) 

5.89 

5.78  (1.11) 

St  at  interp 
(Cd  =  14) 

5.26 

5.17(0.92) 

4 .57 

4  .49  (0.83) 

4.57 

4.50  (0.83) 

Stat  inter? 
(Cd  =  7) 

5.02 

4.95  (0.82) 

4  .28 

4.21  (0.74) 

4.28 

4.21  (0.74) 

Stat  Interp 
(Dspd  Cos) 

4.96 

4.89  (0.83) 

U  22 

4  1 16  (0.75) 

4.22 

4.16(0.75) 

Thin  PI  Spl 
(m  *  2) 

4.92 

4.85  (0.84) 

4.15 

4.09  (0.72) 

4.16 

4.10  (0.72) 

La  pi  So  Spl 
(m  *  3) 

6.08 

5.98  (1.05) 

5.43 

5.33  (1.06) 

5.43 

5.33  (1.06) 

frrke/Grdn 
(3  Pass) 

11.82 

11.55  (2.55) 

11.61 

11 .34  (2.51) 

11.61 

11.34  (2.50) 

pseadoEarnes* 
(2  Pass) 

8.89 

8.70  (  1.  85) 

8  .48 

8.29(1.76) 

8.49 

8.30  (1.77) 

pseudoBarnes' 

(Local) 

7.30 

7.15(1.47) 

6.75 

6  .60  (1.40) 

6.76 

6.62(1.41) 

TABLE  6 


P  *  500,  Theta  =  RANDCP,  Delth  =  BANDOM,  Entries:  FUSE 
rg  3  30,  ro  *  10  Mean 
Nusber  cf  realizations  =  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


Grid-to-cbs: 

Obs-to-grid 

PM  Linear 

Bicub  Spl 

PW  Bicub 

Opt  Interp 

7.00 

6.40 

(Cd  =  1C) 

6.85(1.44) 

6.27  (1.30) 

Local  Ol 

7.48 

6.90 

(Cd  =  10) 

7.33(  1.49) 

6.90  (1.36) 

Barnes* 

9.92 

9.52 

2-Pass 

9.70(2.07) 

9.33  (1.89) 

Barnes* 

8.34 

7.89 

(Local) 

8.21  (1.46) 

7.76  (1.39) 

Stat  Interp 

7.92 

7.41 

(Cd  =  1«) 

7.77  (1.56) 

7.27  (1.43) 

Stat  Interp 

7.56 

7.01 

(Cd  *  7) 

7.41  (1.51) 

6.88  (1.36) 

Stat  Interp 

7.58 

7.04 

(Dtpd  Cos) 

7.43(1.  53) 

6.90  (1.38) 

Thin  PI  Spl 

7.63 

7.06 

(■  =  2) 

Lapl  Sb  spl 

7.47(  1.58) 

6.92  (1.42) 

(n  *  3) 

Frnke/Grdn 
(3  Pass) 

PseudoBarnee ' 
(2  Pass) 

PseudoBarnes* 

(Local) 


TABLE  7 


2  7 


-•  ^  V' 


Analysis 
FUSE  (StDev) 


Bsl  Bicub 


•  v? 


P  =  500,  Iheta  =  RANDOM,  Delth  *  RANDOM,  Entries: 
rg  *  30,  ic  =  10,  rg  (lie)  =  20 
Ruaber  of  realizations  =  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


RMSE  Analysis 
Mean  RMSE (StDev) 


Gr id-to-cbs: 

PH  Linear 

Bicub  Spl 

PH  Bicub 

obs-to-grid 

Opt  Interp 

7.14 

6.47 

(Cd  *  10) 

6.97  (1.52) 

6.32  (1.35) 

Local  01 

7.56 

6.92 

(Cd  *  10) 

7.39(  1.62) 

6.77(1.41) 

Barnes' 

9.64 

9.21 

2- Pass 

9.42(2.04) 

9.01  ( 1 . 9  C ) 

Barnes' 

8.30 

7.73 

(Local) 

7.04  (1.43) 

7.62(1.33) 

Stat  Interp 
(Cd  =  14) 

Stat  Interp 
(Cd  =  7) 

Stat  Interp 
(Dtpd  Ccs) 

Thin  PI  SfI 

7.19 

6.52 

(n  *  2) 

7.04(  1.45) 

6.39  (1.26) 

Bsl  Bicub 


Lapl  Sn  Spl 
(»  *  3) 


Frnke/Grdn 
(3  Pass) 


PseudoBarnes' 
(2  Pass) 

PseudoBarnes' 

(Local) 


TABLE  8 
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E  *  500,  Theta  *  RANDOM,  Delth  *  RANDOM,  Entries:  RMSE  Analysis 
rg  *  2C,  rc  =  10,  rg(lie)  =  30  Mean  RMSE(StDav) 

Rueber  of  realizations  *  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


Gr id-tc-cbs: 

EH  Linear 

Bicub  Spl 

PH  Bicub 

Obs-to-grid 

Opt  Interp 
(Cd  *  1C) 

6.32 

6.17(1.36) 

5.72 

5.60  (1.21) 

Local  oi 
(Cd  =  10) 

6.60 

6.47(  1.33) 

6.05 

5.94  (1.16) 

Barnes* 

2- Eass 

7.10 

6.95(  1.45) 

6.69 

6.55  (1.39) 

Baines* 

(Local) 

7.09 

6.98  (1.22) 

6.54 

6.45  (1. 08) 

Scat  Interp 
(Cd  »  14) 

stat  Interp 
(Cd  *  7) 

Stat  Interp 
(Dipd  Ccc) 

Thin  PI  Spl 

6.31 

5.75 

(m  *  2)  6.17(  1.30)  5.63  (1.16) 

Lapl  Sm  Spl 
(■  =  3) 

Frnke/Grdn 
(3  Eass) 

PsecdoBaines' 

(2  Eass) 

PseudoBaxnes' 

(Local) 


TABLE  9 


E  *  500,  Theta  =  BAN DCM,  Delth  *  RANDOM,  Entries:  BMSE 
rg  *  3C,  rc  =  10,  ro  (lie)  =  5  Mean 
Nuaber  of  realizations  *  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


Gr id-to-cbs: 

EH  Linear 

Obs-to-grid 

Opt  Interp 
(Cd  *  10) 

7.51 

7.37  (1.43) 

Local  01 
(Cd  =  10) 

7.96 

7.85(  1.36) 

Barnes' 

2- Eass 

S  81 

9.60  (  2.  0  1) 

Barnes' 

(Local) 

8.47 

8.36  (1.32) 

stat  Interp 
(Cd  =  14) 

Stat  Interp 
(Cd  «  7) 

Stat  Interp 
(Dipd  Ccs) 

Thin  Pi  Spl 
(■  =  2) 

7.65 

7.51  ( 1.  43) 

Lapl  Sa  Spl 

(■  *  3) 


Bicub  Spl  PH  Bicub 


6.88 

6.74  (1.38) 
7.37 

7.26  (1.27) 
9.28 

9.06  (1.99) 
7.91 

7.80  (1.31) 


6.94 

6.80  (1.39) 


Frnke/Grdn 
(3  Eass) 


PseudoBarnes* 
(2  Eass) 


Analysis 
RMSE (StDev) 


Bsl  Bicub 


PseudoBarnes' 

(Local) 


RUSE  Analysis 
Mean  RMSE(StDev) 


{  -  500,  Theta  =  RANDOM,  Delth  =  RANDOM,  Entries: 
ig  =  30,  ic  =  5 ,  ro(lie)  =10 
Nueber  of  realizations  =  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


Gr id-to-cbs : 

Obs-to-grid 

PH  Linear 

Bicub  Spl 

PH  Bicub 

Opt  Interp 

5.18 

4.19 

(Cd  =  1C) 

5.10(0.91) 

4.13  (0.71) 

Local  Ol 

5.66 

4.76 

(Cd  =  10) 

5.56(1.04) 

4.69  (0.81) 

Barnes' 

9.09 

8.62 

2-Eass 

8.93(  1.69) 

8.47  (1.61) 

Barres' 

6.71 

6.06 

(Local) 

Stat  Interp 
(Cd  =  14) 

Stat  Interp 
(Cd  =  7) 

Stat  Interp 
(Dipd  Ccs) 

6.63  (1.07) 

5.97  (1.04) 

Thin  PI  Spl 

6.13 

5.34 

(n  *  2) 

Lapl  Sm  Spl 
(n  =  3) 

Frnke/Grdn 
(3  lass) 

PseudoBarnes • 
(2  Pass) 

6.02(  1.  15) 

5.25  (0.95) 

PseudoBarnes' 

(Local) 


I 


P  =  50C,  Theta  =  BANDCB,  Delth  =  BANDOM,  Entries: 

rg  =  5 ,  r  c  =  10 

Number  cf  realizations  =  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


BMSE 

Mean 


Grid-to-cbs: 

Obs-to-grid 

PH  Linear 

Bicub  Spl 

PH  Bicub 

Opt  Interp 

3.58 

3.32 

(Cd  =  1C) 

3.44  (1.01) 

3.18  (0.94) 

Local  01 

3.77 

3.53 

(Cd  =  10) 

2.63  (  1.02) 

3.40  (0.  9769 

Barnes* 

4.44 

3.98 

2-Bass 

4.32(  1.01) 

3.87  (0.90) 

Barnes' 

6.26 

5.73 

(Local) 

6.15  (1.16) 

5.64  (1.00) 

Stat  Interp 

3.61 

3.41 

(Cd  *  14) 

3.46  (1.03) 

3.27  (0.96) 

Stat  Interp 

3.65 

3.39 

(C  d  *  7) 

3.50(1.03) 

3.25  (0.  97) 

stat  Interp 

3.78 

3.51 

(Dxpd  Cos) 

3.62(  1.07) 

3.36  (1.02 

Thin  PI  Spl 

4.00 

3.85 

(■  =  2) 

3.84  (  1.  13) 

3.70  (1.07) 

La pi  Sm  Spl 
(■  -  3) 

Frnke/Grdn 
(3  Bass) 

PseudoBarnes* 
(2  Bass) 

PseudoBaxnes1 

(Local) 


Analysis 
BMSE (StDev) 


Bsl  Bicab 


TABLE  12 


E  =  50  0,  Theta  =  RANDOM,  Delth  =  RANDOM,  Entries:  RMSE  Ar.alysi 
ig  =  30,  rc  =  0  Mean  FMSE(St 

Number  of  realizations  =  100 

13X9  grid  of  2.5  degrees,  36  Observation  points 


Grid-to-cts : 

Obs-to-grid 

fw  Linear 

Bicut  Spl 

PH  Bicub 

Opt  Interp 

6.06 

0.74 

(Cd  =  10) 

5.84  (1.  62) 

0.69  (0.27) 

Local  01 

3.86 

2.40 

(Cd  =  10) 

3.79  (0.76) 

2.2. 30(0.69 

Barnes* 

9.38 

8.96 

2-Pass 

9.16  (2.05) 

8.74  (1.98) 

Barnes ' 

6.27 

5.63 

(Local) 

6.  1 5  (  1.21) 

5.49  (1.23) 

St  at  Interp 

7.62 

1.23 

(Cd  =  14) 

7.10  (2.74) 

1  .  16  (0.40) 

Stat  Interp 

3.70 

1.29 

(Cd  =  7) 

3.64  (0.68) 

1.20  (0.  48) 

Stat  Interp 

4.48 

1.03 

(Dmpd  Cos) 

4.37  (0.96) 

0.97  (0.36 

Thin  PI  Spl 

3.74 

2.18 

(m  *  2) 

3.66  (0.78) 

2.02  (0.83) 

Lapl  Sm  Spl 
(m  *  3) 


Frrke/Grdn 
(3  Pass) 

PseudoEarnes' 
(2  Pass) 


Pseudo  Earnes* 
(Local) 


TABLE  13 


in  o 


RMSE  Analysis 
Mean  RMSE(StDrv) 


E  =  500,  Theta  =  100,  Eelth  =  0  Entries: 

rg  *  30,  ro  =  10 

Number  of  realizations  =  100 

11  BY  8  grid  of  5  degrees,  67  Observation  points 


Grid-to-obs : 

Obs-to-grid 

EH  Linear 

Bicub  Spl 

PW  Bicub 

3sl  Bicub 

Opt  Interp 

12.84 

7.62 

7.93 

8.11 

(Cd  =  10) 

12.74  (1.62) 

7.53  (1.19) 

7.84  (1.20) 

8.02(1.19) 

Local  OI 

13.33 

8.44 

8.74 

8.92 

(Cd  =  10) 

13.22  (1.73) 

8.33  (1.  33) 

8.63  (1.35) 

8.82  (1.32) 

Barnes* 

14.33 

10.62 

10.82 

10.95 

2- Pass 

14.21  (1.85) 

10  .49  (1.70) 

10.68  (1.70) 

10.81  (1.71) 

Barnes* 

14.00 

8.82 

9.12 

9.40 

(Local) 

13.91(1.55) 

8.71  (1.36) 

9.79  (1.35) 

9.31  (1.33) 

St  at  interp 

13.75 

8.80 

9.02 

9.20 

(Cd  =  14) 

13.25(  1.57) 

8.70(1.33) 

8.92  (1.31) 

9. 10  (1. 32) 

Stat  Interp 

13.44 

8.31 

8.62 

8.79 

(Cd  =  7) 

13.35  (1.63) 

8.23  (1.18) 

8.53  (1.21) 

8.70(1.19) 

Stat  Interp 

13.57 

8.47 

8.78 

8.95 

(Dmpd  Cos) 

13.47  (1.65) 

8 .38  (1. 25) 

8.68  (1.27) 

8.86  (1.25) 

Thin  PI  Spl 

13.17 

8.08 

8.36 

8.47 

(■  *  2) 

13.07 (1.59) 

7 .99 (1.  17) 

8.27  (1.19) 

8.39  (1. 19) 

Lapl  Sm  Spl 

17.12 

11  .87 

12.08 

12.16 

(m  *  3) 

17.01  (1.86) 

11.76  (1.60) 

11.97  (1.57) 

12.05  (1.58) 

Frnke/Grdn 

17.29 

15.23 

15.31 

15.36 

(3  Eass) 

17.14(2.  29) 

15 .04  (2.35) 

15.13(2.34) 

15.17 (2.37) 

PseudoBarnes 

•14.27 

10.62 

10.82 

10.95 

(2  Eass) 

14.14  (1.87) 

10.49  (1.70) 

10.69  (1.70) 

10.82  (1.71) 

PseudoBarnes 

•13.89 

9.73 

10.00 

10.17 

(Lccal) 

13.77  (1.82) 

9  .60  (1.  58) 

9.87  (1.58) 

10.05  (1.57) 

TABLE  14 


Estimated  and  (empirical)  RMS  errors  for  stati 


method 


Opt  Inter f 
(cd  =  10) 

Local  01 
(cd  =  10) 

Stat  Inteip 
(cd  =  14) 

Stat  InteiF 
(cd  =  7) 

Stat  Interp 
(Dmpd  Cos) 


rg:  5  20  30 

ro:  10  10  10 


3.30 

(3.32) 


5.62  6.29  3 

(5.76)  (6.40)  (3 


3.  55  5.96  6.80 

(3.53)  (6.10)  (6.  90] 


2.90  4.60  5.07 

(3.41)  (6.31)  (7.41] 


3.  76 
(3.39) 

3.82 

(3.51) 


7.14 

(6.25) 

7.00 

(6.32) 


8.26 

(7.011 

7.94 

(7.04] 


TABLE  15 


Figure 


CONTOURS,  THETR  =  100,  DELTH 


Figure 


RMSR/RMS  OBS) 


T3:  RG  -  30,  RO  =  10,  13X9  2.5  DEGREE  GRID 


T5 :  RG  -  20,  R0  -  10,  13X9  2.5  DEGREE  GRID 


30,  R0  -  10,  13X9  2.5  DEGREE  GRID 
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Figure  7 


T14:  RG  -  30,  R0  =  10,  11X8  5  DEGREE  GRID 
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Figure  8 
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PERFORMANCE  DEGRADATION,  RG  -  30,  RO  -  5 
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PERFORMANCE  DEGRADATION,  RG  -  20,  RO  =  10 
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COMMANDER  &  DIRECTOR 
ATTN:  DELAS-D 
U.S.  ARMY  ATMOS.  SCI.  LAB 
WHITE  SAND  MISSILE  RANGE 
WHITE  SANDS,  NM  88002 


DIRECTOR 

NATIONAL  METEORO.  CENTER 
NWS,  NOAA 
WWB  W32,  RM  204 
WASHINGTON,  DC  20233 


CHIEF  OF  NAVAL  RESEARCH  (2) 
LIBRARY  SERVICES,  CODE  784 
BALLSTON  TOWER  #1 
800  QUINCY  ST. 

ARLINGTON,  VA  22217 


COMMANDING  OFFICER 
FLENUMOCEANCEN 
MONTEREY,  CA  93940 


NAVAL  POSTGRADUATE  SCHOOL 
OCEANOGRAPHY  DEPT. 
MONTEREY,  CA  93943 


COMMANDER 

NAVAIRSYSCOM  (AIR-330) 
WASHINGTON,  DC  20361 


SUPERINTENDENT 
ATTN:  USAFA  (DEG) 

USAF  ACADEMY,  CO  80840 


AFGL/OPI 

HANSCOM  AFB,  MA  01731 


DIRECTOR  (2) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA  22314 


ACQUISITIONS  SECT.  IRDB-D823 
LIBRARY  &  INFO.  SERV.,  NOAA 
6009  EXECUTIVE  BLVD. 
ROCKVILLE,  MD  20852 


OFFICE  OF  NAVAL  RESEARCH 
CODE  422AT 
ARLINGTON,  VA  22217 


SUPERINTENDENT 
LIBRARY  REPORTS 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 


LIBRARY  (2) 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


COMMANDER 

NAVOCEANSYSCEN 

DR.  J.  RICHTER,  CODE  532 

SAN  DIEGO,  CA  92152 


AFGWC/DAPL 

OFFUTT  AFB,  NE  68113 


OFFICER  IN  CHARGE 
SERVICE  SCHOOL  COMMAND 
DET.  CHANUTE/STOP  62 
CHANUTE  AFB,  IL  61868 


NOAA-NESDIS  LIAISON 
ATTN:  CODE  SC2 
NASA-JOHNSON  SPACE  CENTER 
HOUSTON,  TX  77058 


DIRECTOR 

OFFICE  OF  PROGRAMS  RX3 
NOAA  RESEARCH  LAB 
BOULDER,  CO  80302 
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NATIONAL  WEATHER  >..aV.  i 
WORLD  WEATHER  B'.!)G. ,  RM  307 
5200  AUTH  ROAD 
CAMP  SPRINGS,  MD  20023 


DIRECTOR 

TECHNIQUES  DEVELOPMENT  LAB 
GRAMAX  BLDG. 

8060  13TH  ST. 

SILVER  SPRING,  MD  20910 


COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN:  DR.  WILLIAM  GRAY 
FORT  COLLINS,  CO  80523 


CHAIRMAN,  METEOROLOGY  DEPT. 
CALIFORNIA  STATE  UNIVERSITY 
SAN  JOSE,  CA  95192 


.)i  •■’LcroR 

NATIONAL  SEVERE  STORMS  LAB 
1313  HALLEY  CIRCLE 
NORMAN,  OK  73069 


LABORATORY  FOR  ATMOS.  SCI. 

NASA  GODDARD  SPACE  FLIGHT  CEN. 
GREENBELT,  MD  20771 


ATMOSPHERIC  SCIENCES  DEPT. 
UCLA 

405  HILGARD  AVE. 

LOS  ANGELES,  CA  90024 


COLORADO  STATE  UNIVERSITY 
ATMOSPHERIC  SCIENCES  DEPT. 
ATTN:  LIBRARIAN 
FT.  COLLINS,  CO  80523 


DIRECTOR 

GEOPHYS.  FLUID  DYNAMICS  LAB 
NOAA,  PRINCETON  UNIVERSITY 
P.O.  BOX  308 
PRINCETON,  NJ  08540 


EXECUTIVE  SECRETARY,  CAO 
SUBCOMMITTEE  ON  ATMOS.  SCI. 
NATIONAL  SCIENCE  FOUNDATION 
RM.  510,  1800  G.  STREET,  NW 
WASHINGTON,  DC  20550 


CHAIRMAN,  METEOROLOGY  DEPT. 
UNIVERSITY  OF  OKLAHOMA 
NORMAN,  OK  73069 


NATIONAL  CENTER  FOR  ATMOS. 

RSCH. ,  LIBRARY  ACQUISITIONS 
P.O.  BOX  3000 
BOULDER,  CO  80302 


UNIVERSITY  OF  WASHINGTON 
ATMOSPHERIC  SCIENCES  DEPT. 
SEATTLE,  WA  98195 


UNIVERSITY  OF  HAWAII 
METEOROLOGY  DEPT. 
2525  CORREA  ROAD 
HONOLULU,  HI  96822 


UNIVERSITY  OF  MARYLAND 
METEOROLOGY  DEPT. 
COLLEGE  PARK,  MD  20742 


CHAIRMAN,  METEOROLOGY  DEPT. 
UNIVERSITY  OF  UTAH 
SALT  LAKE  CITY,  UT  84112 


CHAIRMAN,  METEOROLOGY  DEPT. 
PENNSYLVANIA  STATE  UNIV. 

503  DEIKE  BLDG. 

UNIVERSITY  PARK,  PA  16802 


DIRECTOR 

COASTAL  STUDIES  INSTITUTE 
LOUISIANA  STATE  UNIVERSITY 
ATTN:  0.  HUH 
BATON  ROUGE,  LA  70803 


CHAIRMAN 

ATMOS.  SCIENCES  DEPT. 
UNIVERSITY  OF  VIRGINIA 
CHARLOTTESVILLE,  VA  22903 


ATMOSPHERIC  SCIENCES  CENTER 
DESERT  RESEARCH  INSTITUTE 
P.O.  BOX  60220 
RENO,  NV  89506 
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FLORIDA  STATE  UNIVERSITY 
ENVIRONMENTAL  SCIENCES  DEPT. 
TALLAHASSEE,  FL  32306 


ATMOSPHERIC  SCIENCES  DEPT. 
OREGON  STATE  UNIVERSITY 
CORVALLIS,  OR  97331 


CHAIRMAN 

METEOROLOGY  DEPT. 
MASSACHUSETTS  INSTITUTE  OF 
TECHNOLOGY 
CAMBRIDGE,  MA  02139 


ATMOSPHERIC  SCI.  RSCH.  CENTER 
NEW  YORK  STATE  UNIV. 

1400  WASHINGTON  AVE. 

ALBANY,  NY  12222 


SYSTEMS  &  APPLIED  SCI.  CORP.  CENTER  FOR  ENV.  &  MAN,  INC. 
ATTN:  LIBRARY,  SUITE  500  RESEARCH  LIBRARY 

6811  KENILWORTH  AVE.  275  WINDSOR  ST. 

RIVERDALE,  MD  20840  HARTFORD,  CT  06120 


METEOROLOGY  INTL.,  INC. 
P.O.  BOX  22920 
CARMEL,  CA  93922 


OCEAN  DATA  Sv STEMS,  INC. 
2460  GARDEN  ROAD 
MONTEREY,  CA  93940 


DEPT.  OF  MATHEMATICS 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


DR.  RICHARD  LAU 
OFFICE  OF  NAVAL  RESEARCH 
1030  EAST  GREEN  ST. 
PASADENA,  CA  91106 


LIBRARY,  FLEET  NUMERICAL( 2) 
OCEANOGRAPHY  CENTER 
MONTEREY,  CA  93943 


PROFESSOR  W.J.  GORDON 
CENTER  FOR  SCIENTIFIC 
COMPUTATION  &  INTERACTIVE 
GRAPHICS 

DREXEL  UNIVERSITY 
PHILADELPHIA,  PA  19104 


PROFESSOR  DONALD  GAVER 
DEPARTMENT  OF  OPERATIONS 
RESEARCH 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


CONTROL  DATA  CORP. 
METEOROLOGY  DEPT.  RSCH.  DIV. 
2800  E.  OLD  SHAKOPEE  RD. 

BOX  1249 

MINNEAPOLIS,  MN  55440 


EUROPEAN  CENTRE  FOR  MEDIUM 
RANGE  WEATHER  FORECASTS 
SHINFIELD  PARK,  READING 
BERKSHIRE  RG29AX,  ENGLAND 


PROFESSOR  G.  LATTA,  CODE  53LZ 
CHRMN,  DEPARTMENT  OF  MATH 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


PROFESSOR  R.E.  BARNHILL 
DEPARTMENT  OF  MATHEMATICS 
UNIVERSITY  OF  UTAH 
SALT  LAKE  CITY,  UT  84112 


DR.  EDWARD  BARKrr>  (10) 

NAVAL  ENVIRONMEK...L.  PREDICTION 
RESEARCH  FACILITY 
MONTEREY,  CA  93943 


MS.  LINDA  THIEL 
DEPARTMENT  OF  MATHEMATICAL 
SCIENCES 

DREXEL  UNIVERSITY 
PHILADELPHIA,  PA  19104 


PROFESSOR  P.A.  JACOBS 
DEPARTMENT  OF  OPERATIONS 
RESEARCH 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


METEOROLOGY  RESEARCH,  INC. 
464  W.  WOODBURY  RD. 
ALTADENA,  CA  91001 


SCIENCE  APPLICATIONS,  INC. 

205  MONTECITO  AVENUE 

MONTEREY,  CA  93940 


DEAN  OF  RESEARCH  (2) 
CODE  012 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


PROFESSOR  R.  FRANKE,  CODE  53FE( 1 5 
DEPARTMENT  OF  MATHEMATICS 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


PROFESSOR  G.M.  NIELSON 
DEPARTMENT  OF  MATHEMATICS 
ARIZONA  STATE  UNIVERSITY 
TEMPE ,  ARIZONA  85281 


DR.  THOMAS  ROSMOND 
NAVAL  ENVIRONMENTAL  PREDICTION 
RESEARCH  FACILITY 
MONTEREY,  CA  93943 


PROFESSOR  GR."CE  WAHBA 
DEPARTMENT  OF  STATISTICS 
UNIVERSITY  OF  WISCONSIN 
MADISON,  WI  53705 


PROFESSOR  H.  JEAN  THIEBAUX 
DEPARTMENT  OF  MATHEMATICS, 
STATISTICS  &  COMPUTING  SCIENCE 
DALHAUSIE  UNIVERSITY 
HALIFAX,  NOVA  SCOTIA, 

CANADA  133H  4H8 


OR.  R.S.  SEAMAN 
AUSTRALIAN  NUMERICAL 
METEOROLOGY  RESEARCH  CENTRE 
P.O.  BOX  5089AA 
MELBOURNE,  VICTORIA, 
AUSTRALIA,  3001 


DR.  LAWRENCE  BREAKER 
DEPARTMENT  OF  OCEANOGRAPHY 
NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943 


PROFESSOR  PETER  ALFELD 
DEPARTMENT  OF  MATHEMATICS 
UNIVERSITY  OF  UTAH 
SALT  LAKE  CITY,  UT  84112 


PROFESSOR  MARK  E.  HAWLEY 
DEPARTMENT  OF  ENVIRONMENTAL 
SCIENCES 

UNIVERSITY  OF  VIRGINIA 
CHARLOTTESVILLE,  VA  22903 


